Dynamic saturation in semiconductor 
optical amplifiers: accurate model, role 
of carrier density, and slow light 



Perrine Berger 1 2 , Mehdi Alouini 1 3 , Jerome Bourderionnet 1 , 
Fabien Bretenaker 2 , and Daniel Dolfl 1 

1 Thales Research & Technology, I av. Augustin Fresnel, 91 767 Pcdaiseau Cedex, France 

2 Laboratoire Aime Cotton, CNRS-Universite Paris Sud 11, Campus d'Orsay, 91405 Orsay 

Cedex, France 

3 Institut de Physique de Rennes, UMR CNRS 6251, Campus de Beaulieu, 35042 Rennes 

Cedex, France 
Perrine.Berger@thalesgroup.com 



Abstract: We developed an improved model in order to predict the RF 
behavior and the slow light properties of the SOA valid for any experimental 
conditions. It takes into account the dynamic saturation of the SOA, which 
can be fully characterized by a simple measurement, and only relies 
on material fitting parameters, independent of the optical intensity and 
the injected current. The present model is validated by showing a good 
agreement with experiments for small and large modulation indices. 
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1. Introduction 

The generation of continuously tunable optical delays is a key element in microwave photon- 
ics. Among the targeted applications, one can quote the filtering of microwave signals, the 
synchronization of optoelectronics oscillators, and the control of optically fed phased array an- 
tennas Ii||2][3). With these applications in view, large efforts are currently done in order to 
develop delay lines based on slow and fast light effects J4j |5J 13 HI • To date, one of the most 
mature approaches for integration in real field systems is that based on Coherent Population 
Oscillations (CPO) in semiconductor structures ll9l ITOl fTTTl . This approach offers compactness, 
continuous tunability of the delay through injected current control, and possible high-level par- 
allelism lfT2l[T3l . Obviously, the implementation of CPO effects in microwave photonics delay 
lines relies on accurate theoretical description of the underlying mechanisms in order to develop 
reliable predictive models. Numerous theoretical models have been developed in the past few 
years to describe CPO effects in Semiconductor Optical Amplifiers (SOAs) (14\ [TSl [161 [T7l . 



They are usually based on a semi-classical description of the interaction between the carriers 
and the input optical fields. These models offer a comprehensive understanding of the gain sat- 
uration dynamics and associated group index changes. However, on the one hand, a complete 
model would require a detailed knowledge of the geometrical and material parameters of the 
semiconductor structure lfl8l \19\ . Unfortunately, most of them are unknown especially when 
the SOA under consideration is a commercially available device. On the other hand, others, 
simpler, assumed that both the saturation power and the carrier recombination lifetime are con- 
stant (5] [20] Ell ■ This assumption applies when the SOA is operated at a fixed injection 
current j22l 23|. However, the injection current and the input optical power have to be tuned 
over a wide range in order to control the speed of light into the SOA: consequently this assump- 
tion restricts the predictive capability of a model describing microwave-photonics delay lines 
using slow light in SOA. 

In this paper we derive an improved model that enables to predict the RF gain compression, 
the RF phase delay, and the optical group delay and which is valid for all experimental condi- 
tions for a given component. Furthermore, we show that the detailed knowledge of the inner 
geometrical and material characteristics of the SOA is not required provided that some prelimi- 
nary and easy characterization measurements are conducted. This model is then experimentally 
validated. 



2. Model 

We consider an optical carrier modulated by an RF signal and injected in a traveling wave SOA. 
The total field is then composed of the optical carrier of complex amplitude Eq and two side- 
bands of complex amplitudes E\ and E%. The total optical field E is normalized to include the 
factor y/EQiiQCQ, i.e., the optical intensity is given by I opt (z,t) = k\Etotai\ 2 = U +Me~' n ' +c.c, 
under small RF signal approximation. U is the DC component of the intensity, M = 
\ (EqE% +E\Eq) is the beat-note term at the RF frequency ii. 

The local equations for the propagation of the optical field E totc ,i and the evolution of carrier 
density N inside the SOA are ifPfll : 



dN(z,t) _ I N(z,t) g(z,t)\E total (z,t)l 



dt qV T v hCO 

\E tota ,(z,t)\ 2 [- 7 +r g (z,t)}, (2) 



d\E lula ,(z,t)\ 2 



dz 

where 7 holds for the internal losses of the SOA, Fg(z,t) is the material modal gain, t s is 
the carrier lifetime, / is the injected current, V is the volume of the active region, and (0 is 
the pulsation of the optical carrier Eq. We introduce N(z,t) — N(z) + AN(z)e~' nt +c.c. and 
g(z,t) = g(N(z,t)) = g(N(z))+a(N(z))AN(z,t)e- iQt + c.c. where a is the differential gain 
a(N) = ^||;v. The wavelength of the optical carrier is fixed. Consequently, the equations Eq.Q] 
and Eq.|2]lead to: 



^ = U[- Y +Tg(N)], (3) 

^ = M(- y+ rWi u { u jW n (4) 

dz \ * SV >\ l + U/U s (N)-ito s (N)J y 

where U S (N) is the saturation intensity defined as: U S (N) = Q ^'"^ ■ 

In most of the simple models, the common approach to solve equations d3} and is to 
consider a and T s constant with respect to the carrier density and thus over the whole length 



of the device |5] [14] |20l |2T| . This approximation does however not give account of strong 
saturation conditions, with high gain and carrier density variations, which typically occur in 
quantum wells structures with strong carrier confinement. In this paper, we propose to consider 
the carrier density variation along the propagation axis and its influence on a and T s , Our central 
hypothesis is that a and T s can be determined as functions of the DC component of the optical 
intensity U solely, allowing these dependencies to be determined from gain measurements. 

Let us first suppose that we fulfill the small signal condition. In this case, the stimulated 
emission is negligible compared to the spontaneous emission, leading to the unsaturated steady 
state solution of the rate equation (Eq.[T|): 

' s (5) 



Q L S ac t T s 

where L is the length of the SOA, S aa is the area of the active section of the SOA. Moreover, 
we also suppose in this case that the carrier density N is constant along the SOA. These hy- 
pothesis are equivalent to consider that the amplified spontaneous emission does not saturate 
the gain. A verification of this assumption will be shown in section[3] Under these conditions, a 
measurement of the small signal modal gain Fgo versus / will be equivalent, owing to Eq.[5] to 
a determination of the modal gain Fg versus N / T s . Here, T is the ratio S act /S gU id e of the active 
to modal gain areas in the SOA. 

A last relationship between ^ and U is then required to determine the modal gain Fg as a 

function of U. It is obtained by substituting Fg(^-) in the saturated steady state solution of the 
carriers rate equation (Eq.Q~|i: 

' * 0, (6, 



q L Sact Tj fl(0 F 

where the injected current / is now fixed by the operating conditions. 

Added to the previous relationship between Fg and j-, the Eq.|6]gives another expression of 
Fg as a function of ^, ^ and /. Consequently, Fg and j. can be known with respect to the local 
intensity and the injected current /. 

To solve Eq.|31 we need to express N as a function of ^4 and /. This is equivalent to express 

with respect to j. since j. is known as a function of and /. Consequently, we model our 
SOA using the well-known equation 



— =AN + BN 2 + CN 3 , (7) 

where A, B, and C, which are respectively the non-radiative, spontaneous and Auger recom- 
bination coefficients, are the only parameters that will have to be fitted from the experimental 
results. 

Using Eq.|7]and the fact that we have proved that N/x s and Fg can be considered as function 
of an d I on ly> we see ma t N, Fa = F^, and ^ = can also be considered as functions 
of an d I- This permits to replace Eqs. (O and (HJi by the following system: 
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Fig. 1. (a) Experimental fiber-to-fiber gain G with respect to the output optical power P ou , 
at a strong current (500 mA). The double arrow indicates the range of the measured ASE 
output power, (b) Experimental small signal gain Tgo as a function of the injected current / 
at 1535 nm, and fitted by: Tg = C\ - % with C\ = 5588.7m" 1 and C 2 = 306. lA^.m" 1 . 
(c) Deduced material modal gain Tg(U) as a function of the local intensity U at 500 mA. 



Eqs. (O and (O are then numerically solved: Eq. [8] gives with the initial condi- 

tion = v^FT' w here P,„ is the optical input power. can be then introduced into 
Eq. |9j and the microwave transfer function of the SOA S21 = Tf^^y. where j\ are the in- 
sertion losses, is then computed. If the output power of the RF microwave signal is wanted, 
Puf = 2Rrip h Yj\ mP >"M( L ) w i tn tne m i);i a i condition M(0) = 1, where m the input modulation 
index, and R and rjp/, are respectively the photodiode resistive load and efficiency. 

The microwave complex transfer function S21 fully characterizes the slow light properties of 
the SOA. Indeed the optical group delay Ax g can be expressed as Az g (Cl) = " , and the 

group index An g {£l) = - 

It is important to note that the recombination coefficients A, B and C are the only fitting 
parameters of our model. Once obtained from experimental data, they are fixed for any other 
experimental conditions. Moreover, the only geometrical required parameters are the length 
L of the SOA and the active area cross section S uc t- The derivation of a predictive model, 
independent of the experimental conditions (current and input optical power) is then possible, 
provided that the simple measurements of the total losses and the small signal gain versus 
the current are conducted. The above model lies in the fact that first, the spatial variations of 
the saturation parameters are taken into account, and second, their values with respect to the 
local optical power are deduced from a simple measurement. These keys ideas lead to a very 
convenient model of the microwave complex transfer function of the SOA, and then of the slow 
light properties of the component. It can be easily used to characterize commercial components 
whose design details are usually unknown, as we will experimentally show in the next section. 

3. Experiment 

In order to validate our model, we studied a commercially available SOA (InP/InGaAsP Quan- 
tum Well Booster Amplifier from COVEGA). The length L of this SOA is 1.50 mm and the 
active area cross-section is set at 0.06 )J.m 2 . We proposed to compare the experimental and 
simulated complex transfer function S21 for a large set of operating conditions (Pj n ,I). As ex- 
plained in section|2j the study of the phase arg(S2\ ) is equivalent to the optical group delay. We 
can then restrain our comparison to S21 ■ In order to fully characterize the response of the SOA 
through Tg{U) as described in section the preliminary step consists in measuring the total 
losses and the unsaturated gain Egot/) for different injected currents. 

The total losses are measured by the following experiment: at low current, the output optical 
power is measured while a strong input optical power is sent into the SOA. When the current 
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Fig. 2. Experimental set-up. For small modulation index m, a laser is externally modulated 
by a Mach-Zehnder modulator (MZ) (a); for large modulation index, a directly modulated 
laser is used (b). In both cases, the input optical power P,„ is controlled through a variable 
optical attenuator; two optical isolators are used before and after the SOA. The photode- 
tector (PD) restitute the RF signal. The Vector Network Analyser (VNA) is calibrated with 
the whole link without the SOA, in order to measure the RF transfer function of the SOA. 



is low enough, the SOA is in the absorption regime: the resulting output power is an increasing 
function of the input power (absorption saturation). When the current is above the transparency 
current, the resulting output power becomes a decreasing function of the input power (gain 
saturation). Between these two regimes, i.e. at transparency, the ratio between the output power 
and the input power is exactly equal to the total losses. The total losses of the SOA ^; exp(— yL) 
are measured to be equal to — 16.4dB in our case. 

To measure only the unsaturated gain Fgo despite the amplified spontaneous emission, we 
measured the SOA unsaturated RF gain Grf (=1^21 1 2 ) at a RF frequency Q. well above 1/t v 
(typically 20 GHz). The derivation of the modal gain Fg with respect to j- from the unsaturated 
gain Tgo(I) is relying on the hypothesis that the amplified stimulated emission (ASE) does not 
saturate the gain. In Fig. we represent the experimental fiber-to-fiber gain with respect to 
the output optical power at a strong current (500 mA) and the range of the experimental output 
power of the ASE. The maximum power of the ASE is equal to 1 .54 dBm. Moreover, when the 
small signal measurement is performed, a maximum input optical power of 80 jiW was used, 
corresponding to an output optical power of 8.1 dBm for the maximum current. Consequently, 
both signal and ASE output power level are well below the output power required to saturate 
the gain (14.2 dBm for a 3dB gain reduction). Therefore, the experimental conditions match 
our preliminary assumptions. Under these conditions, Eq. [8] can be simplified and integrated, 
leading to the expression of the optical small signal gain Tgo'- exp(TgoL) = - j^Pjj^q ■ As shown 
in Fig. [2^, we used a Vector Network Analyzer (VNA) to measure the RF gain Grf for a small 
input power which does not saturate the SOA (typically 10 — 80/iW). The unsaturated gain of 
our SOA is displayed in Fig.QJ). It is empirically fitted by Tgo = Ci — with a good agreement. 
From this simple measurement and using Eq. [6] the material modal gain Yg is then known as a 
function of the local intensity U /T inside the SOA (Fig.QJ). 

The complex RF transfer function of the SOA is measured thanks to a VNA for small and 
large modulation indices (set-ups in Fig. [2j. In Fig. [3] and Fig.|H we report the corresponding 
RF gains, 20 log \Sji |, and the measured evolution of the RF phase shift, arg(52i ), as a function 
of the modulation frequency i2. In each of these figures, the plots labeled (a) and (b) correspond 
to the evolutions of the RF gains and phase shifts versus RF frequency, for different injected 
currents, while the plots labeled (c) and (d) are obtained by managing the input optical power. 
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Fig. 3. Low modulation index (m = 0.06) : gain and phase shift simulations (dashed line) 
and experimental data (solid line) for (a) and (b): different injected currents at P,„ = OdBm, 
and for (c) and (d): different optical input powers P,,, at / = 500mA. The operating wave- 
length was 1535nm. 

4. Discussion 

In Fig. [3] and Fig. [4] the simulation results are reported in dashed line. The best 
fit values for the recombination coefficients are: A = 2x 10 9 s~', B~1.2x 10 _10 cm 3 s _1 , 
C = 1.8 x 10 cm 6 s -1 . These values are in the range of what can be found in the literature 
for semiconductor materials l25l l26l |27] [28l . The computed complex transfer function shows 
a very good agreement with the experimental data, both at small and large modulation index, 
for any experimental conditions (injected current, input optical power), and with a single set 
of the fitting parameters (A, B, C): our convenient model is predictive for any experimental 
conditions. 

In order to highlight the weight of the spatial variations of the carrier density and the satura- 
tion parameters, we plotted in Fig . the variations of the carrier density N along the SOA 
for the different experimental situations of Fig. [3] The subsequent variations of the modal gain 
Fg and the saturation parameters P s , T s and a, with respect to N, are displayed in Fig. |5};l5]l 
We find at least one order of magnitude of variation for almost all these parameters, which are 
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Fig. 4. Large modulation index (m > 0.6) : gain and phase shift simulations (dashed line) 
and experimental data (solid line) for (a) and (b): different injected currents at P ln = OdBm. 
The operating wavelength was 1548. 5nm. 

nevertheless often taken constant in literature for practical models 151 [141 l20l |2TI . According 
to Eq. |71 this approximation can be justified when the variations of N along the SOA are rela- 
tively not too strong, that is for moderate bias current (< 150 mA in our case) or a high bias 
current, but low optical power. However, for any other condition, and especially in the case of 
quantum well or quantum dots structures, it is necessary to take into account the saturation dy- 
namics along the propagation to ensure good performances of the model and robustness versus 
changes in experimental conditions. Indeed, Fig.[5jshows that considering P s , z s and a constant, 
and then Yg linear with N, drastically limits the range of experimental conditions (P;„ , /) where 
such models are valid, which forces the saturation parameters to be adjusted with the current 
and/or the optical input power. 

Our improved model is still easy to use, even for commercial components, but despite the 
hypothesis we were compelled to make, it remains valid for a large range of experimental 
conditions, with a reduced set of unknown - and thus fitted- parameters. These advantages have 
been achieved by taking into account the spatial variation of the saturation parameters and by 
showing that their values as a function of the local optical power can be retrieved from a simple 
measurement. It ensures that the model only relies on material fitting parameters, independent 
of the optical intensity and injected current. 

The slow light properties are then also modeled for a large range of the input optical powers 
Pi„ and injected currents /, which is essential from the operational point of view, since the speed 
of light in SOA is controlled by these two key parameters. While the applications of slow light 
in SOA are taking shape, a convenient and accurate model with the parameters tuning the delays 
is a necessary tool to fully characterize the effect of slow light in SOA on a microwave link, or 
to develop new architectures improving the slow light properties. This model could be easily 
used when an optical filtering is performed after the SOA to enhance the slow light effect, as 
described in [20., 29). In this case, Eq. [9]just has to be replaced by the corresponding coupled 
equations in E\ and Ei. Moreover, to take into account higher order coherent population os- 
cillations ll30l . the present model can be generalized using equations similar to Eq.|9]for each 
harmonic of the optical intensity. The determination of Fg as a function of U is slightly more 
subtle in this case: it is presented in another paper, in order to study the harmonic generation 
and the intermodulation products |3D . 
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Fig. 5. (a) and (b): Simulated carrier density N along the SOA: (a) at a fixed input optical 
power (0 dBm), for various currents; (b) at a fixed current (500 mA), for various input 
optical power, (c) and (d): Simulated variations with respect to the carrier density N of 
(c) the modal gain Tg (solid line), and the modal differential gain a (dashed line); (d) the 
carrier lifetime t s (solid line), and the local saturation power P s (dashed line). 



5. Conclusion 

We developed an improved but still convenient model in order to predict the RF behavior and 
slow light properties of the SOA, valid for any experimental conditions (input optical power, 
injected current). It takes into account the spatial variations of the saturation parameters along 
the SOA, which are fully characterized by the simple measurement of the small signal gain. The 
resulting model only relies on material fitting parameters, independent of the optical intensity 
and injected current. We showed a remarkably good agreement between the model and the 
experimental data, at small and large modulation indices. The ease of use and the accurate 
prediction obtained for any experimental conditions will be useful to characterize the effect 
of slow light in SOA on a microwave link, and to develop new designs improving the slow 
light properties. The key ideas of this improved model can easily be used when optical filtering 
is performed after the SOA. A generalization of our approach will be carried out in a next 
step, in order to determine the harmonic generation, intermodulation products and spurious 
free dynamic range, for a full characterization of a SOA based opto-electronic link. 
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